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Interchange Method in Compressible Magnetized Couette Flow: 
Magnetorotational and Magnetoconvective Instabilities 
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ABSTRACT 

We obtain the general forms of the axisymmetric stability criteria in a mag¬ 
netized compressible Couette flow using an energy variational principle, the so- 
called interchange or Chandrasekhar’s method, which we applied successfully in 
the incompressible case. This formulation accounts for the simultaneous presence 
of gravity, rotation, a toroidal magnetic held, a weak axial magnetic held, entropy 
gradients, and density gradients in the initial equilibrium state. The power of 
the method lies in its simplicity which allows us to derive extremely compact and 
physically clear expressions for the relevant stability criteria despite the inclusion 
of so many physical ehects. In the implementation of the method, all the appli¬ 
cable conservation laws are explicitly taken into account during the variations of 
a quantity with dimensions of energy which we call the “free-energy function.” 

As in the incompressible case, the presence of an axial held invalidates the 
conservation laws of angular momentum and azimuthal magnetic hux and intro¬ 
duces instead isorotation and axial current conservation along held lines. Our 
results are therefore markedly diherent depending on whether an axial magnetic 
held is present, and generalize in two simple expressions all previously known, 
partial stability criteria for the appearance of magnetorotational instability. Fur¬ 
thermore, the coupling between magnetic tension and buoyancy and its inhuence 
to the dynamics of nonhomoentropic magnetized hows becomes quite clear from 
our results. In the limits of plane-parallel atmospheres and homoentropic hows, 
our formulation easily recovers the stability criteria for suppression of convective 
and Parker instabilities, as well as some related special cases studied over 40 
years ago by Newcomb and Tserkovnikov via laborious variational techniques. 
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1. Introduction 
1.1. Motivation 

There has been a recent controversy about the direction of transport of angular momen- 
tumin in radiatively inefficient accretion disk theory, the theory most successful in address¬ 
ing the dynamics and radiation emission in low-luminosity accretion flows. Two schools of 
thought have come to different conclusions by analyzing two opposite limiting cases: Pro¬ 
ponents of magnetorotational instability (Balbus & Hawley 1998, 2002; Hawley et ah 2001) 
largely ignore convection and hnd generic outward transport of angular momentum, while 
proponents of convective instability largely ignore magnetic helds and hnd on average slow 
inward transport of angular momentum by the combination of viscous forces and buoyancy 
(Igumenshchev & Abramowicz 1999, 2000; Stone et ah 1999; Igumenshchev et ah 2000; 
Narayan et ah 2000; Quataert & Gruzinov 2000; Abramowicz et ah 2002). 

It is our opinion that this controversy stems from our ignorance of the precise stability 
criterion of nonhomoentropic, magnetized, differentially rotating hows that are thought to be 
associated with the aforementioned low-efficiency accretion onto black holes. Such a criterion 
has eluded standard linearization and energy variational techniques because the correspond¬ 
ing set of perturbed axisymmetric equations is quite complex and the exact problem has 
never been reduced to a single modal equation. For this reason and in view of its current im¬ 
portance to the understanding of low-luminosity accretion onto black holes in X-ray binaries 
and active galactic nuclei (see Narayan [2002] for a review), we revisit this problem employ¬ 
ing the powerful free-energy variational method developed by Christodoulou et al. (1995a,b) 
and applied to incompressible Couette hows in an earlier publication (Christodoulou et al. 
1996, hereafter Paper I). 

Very recently, Balbus & Hawley (2002) and Narayan et al. (2002) attempted to tackle 
the same problem with a set of linearized equations derived by Balbus & Hawley (1991) by 
augmenting them with an additional radial convective term. However, it has been pointed 
out by Knobloch (1992) and re-iterated in Paper I that the original set of equations does 
not take into account the ehect of a curvature-dependent magnetic-tension term. But it is 
precisely this term that couples into the conventional Brunt-Vaisala frequency and modihes 
the ehect of buoyancy in magnetized nonhomoentropic hows. As we shall see below, the 
inhuence of the magnetic tension is particularly dramatic in hows with combined toroidal 
and poloidal magnetic helds, where it can single-handedly change the sign of the convective 
term in the stability criterion, thereby reversing the ehect of superadiabatic and subadiabatic 
entropy gradients relative to what is known from the conventional Schwarzschild criterion. 

While the results of our study shed some light to the above-mentioned controversy, we 
believe that they are also of considerably broader interest. We expect the results described 
in this paper to prove useful in undestanding both the global structure of such accretion 
hows and the associated jet-like outhows; as well as in other areas of astrophysics and likely 
in branches of plasma/fusion physics which involve the combined ehects of convection, diher- 
ential rotation, and magnetic helds. This generality is not surprising when considering that 
some of the hydromagnetic stability criteria discussed below have been previously derived in 
studies of laboratory plasmas and planetary magnetospheres (see, e.g., Newcomb [1962] and 
Rogers & Sonnerup [1986], respectively). 



1.2. Previous Work 


In Paper I we analyzed the stability of magnetized incompressible Couette flows using 
an energy variational principle, the so-called interchange or Chandrasekhar’s (1960) method. 
This formulation reproduced previously known stability criteria in such an elegant and sim¬ 
ple manner, that we became convinced that the interchange method is indeed physically 
correct, mathematically robust, efficient, and an ideal technique for analyzing a variety of 
difficult stability problems. Motivated by this success we proceed in this work to extend the 
interchange method to the stability of axisymmetric compressible nonhomoentropic magne¬ 
tized Couette flows. In doing so, it becomes necessary to clarify a number of crucial details 
concerning the application of the method: Why do the ’’isotropic” pressure forces not con¬ 
tribute to the construction of the free-energy function? What is the role of gravity during 
the simultaneous interchange of two fluid elements? What are the rules according to which 
the interchange of compressible fluid elements is performed? Does the interchange involve 
elements of equal mass or equal volume? 

Understanding the answers to these questions sheds light to the physics of the studied 
instabilities as well as to the physical meaning of the interchange method itself. For exam¬ 
ple, it becomes clear why we should simultaneously interchange two fluid elements and not 
simply displace/perturb one element from its original position. Such a single-displacement 
approach was discussed by Rogers & Sonnerup (1986) for the case of uniform rotation in an 
homoentropic fluid and led to results which did not allow one to discern the basic physics 
of the stability problem. Rogers & Sonnerup also used the standard interchange method of 
two fluid elements, but they over-constrained the volumes of the fluid elements attempting 
to interchange equal masses. Their result is valid for uniform rotation, but it is doubtful 
whether the same methodology can be productive in more complex (differentially rotating 
and nonhomoentropic) fluids. 

Interestingly enough, Rogers & Sonnerup (1986) were trying to understand why a cal¬ 
culation of the stability of magnetized compressible flows previously performed by Cheng 
(1985) had led to two different stability criteria depending on the use of a purely convective 
stability criterion or the interchange method. Cheng argued that the interchange method 
produces incorrect results. Unfortunately, the disagreement is entirely due to the incorrect 
implementation of the interchange method in Cheng’s investigation. In particular, Cheng 
demanded that the interchanged fluid elements be in pressure balance after the interchange. 
This condition is not appropriate which explains why the interchange method produced 
the wrong result. The inconsistency in assuming pressure balance after the interchange is 
easily identified in the “free-energy” variational framework defined by Christodoulou et al. 
(1995a,b): When fluid elements are interchanged, the system is effectively taken out of equi¬ 
librium and the only relevant question is whether the free energy increases or decreases. As 
was graphically and analytically shown in Christodoulou et al. (1995a), a decrease in the 
out-of-equilibrium free energy indicates the presence of another equilibrium state of lower 
energy and the system will evolve “downhill” in energy space seeking this new equilibrium 
state. Whether the system will get there on a dynamical timescale (or in an Alfven time 
as Cheng remarked) depends purely on the applicable conservation laws, but this is of no 
consequence to the variational method and certainly does not imply that the anticipated 
pressure balance in the final equilibrium state should be built into the interchange method. 

Newcomb (1962) has also analyzed the same problem in the special case of a non¬ 
gravitating homoentropic fluid by using a Hamiltonian formulation along with an energy 
variational principle and was able to recognize the importance of the conservation law of 
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angular momentum in the case of a Couette flow threaded only by a toroidal magnetic held 
and the importance of isorotation in the case of a Couette how threaded only by a weak axial 
held. Unfortunately, Newcomb’s method is rather complicated and mathematically tedious, 
a fact that quickly obscures the physics of the problem. These difficulties have prevented 
Newcomb (1962) from obtaining a stability criterion in the combined case of both toroidal 
and weak axial magnetic-held components. We discuss in more detail the above-mentioned 
older investigations in an Appendix to this paper, where we analyze various special cases of 
our general results. 


1.3. Outline 

In the following sections, we describe the results that we have obtained on the above sta¬ 
bility problems by our free-energy variational method and we clarify the role of various force 
components and conservation laws in the general case of rotating, gravitating, magnetized, 
compressible, nonhomoentropic hows. In § 2, we discuss details of the implementation of the 
interchange method to the axisymmetric stability of compressible hows. In § 3, we analyze 
the case of Couette how threaded by a purely toroidal magnetic held. In § 4, we demonstrate 
that the introduction of a weak axial-held component to the model of § 3 ehects formidable 
changes to the applicable conservation laws and to the resulting stability characteristics of 
the how. In § 5, we apply our method to static plane-parallel atmospheres under the inhu- 
ence of gravity and we recover well-known criteria (Schwarzschild 1906; Tserkovnikov 1960; 
Newcomb 1961; Parker 1966, 1975, 1979) for suppression of convective motions and magnetic 
buoyancy. Finally, in § 6, we discuss astrophysical applications of our results in the context 
of accretion hows and magnetically-driven jet outhows. In an Appendix to the paper, we 
reduce our results by selectively ignoring some of the physics involved (entropy gradients and 
either gravity or diherential rotation), and we recover all previously known partial stability 
criteria (Tserkovnikov 1960; Newcomb 1962; Rogers & Sonnerup 1986). 


2. The Interchange Method 

Let us hrst describe the interchange method for compressible hows. We adopt cylindri¬ 
cal coordinates (r, 0, z) and consider axisymmetric perturbations applied to axisymmetric 
((9/(90 = 0) magnetized Couette hows. The magnetic held is assumed to be frozen-in to the 
huid. Following Rayleigh (see, e.g., Chandrasekhar 1981, § 66), we consider the simultane¬ 
ous radial interchange of two nearby huid elements orbiting initially at radii ri and r 2 > ri. 
In order to perform the interchange correctly in the general compressible case, we need to 
account for the following important details: 

1. The initial conhguration is in hydro magnetic equilibrium. 

2. The interchange involves two huid elements of equal volume, not equal mass. During 
and at the end of the interchange, the ambient huid must retain its continuity and 
the method should not perturb or modify any boundaries that may exist in the huid. 
Therefore, we should make sure that each displaced huid element occupies exactly the 
volume left vacant by the displacement of the other huid element. 
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3. Volume continuity also implies that, although the sum of the two masses of the inter¬ 
changed elements is conserved, individual masses are not. Therefore, unequal masses 
are interchanged which is also necessary in the presence of gravity, if the fluid elements 
are going to feel the effects of the gravitational field. 

4. The interchange is performed adiabatically, i.e., in a way that certain physical pa¬ 
rameters which express conservation laws associated with each fluid element remain 
constant during the interchange. Where applicable, the conservation laws must be 
expressed per unit mass since individual element masses are not conserved. 

5. The free-energy variational principle is applied between the initial equilibrium state 
and the new nonequilibrium end state. What will happen afterwards, e.g. whether 
the interchanged fluid elements will expand or contract in their new positions and will 
disturb the ambient fluid locally, is of no consequence to the method. 


It is physically important to realize that if the change in free energy after the interchage 
is positive (negative), then the initial equilibrium is dynamically stable (unstable) to axisym- 
metric perturbations and the exchanged fluid elements are subject to return to (escape from) 
their original equilibrium positions (cf. Christodoulou et al. [1995a,b]). In implementing the 
method, we hrst calculate the initial condition for equilibrium and the free energy spent or 
gained due to the action of radial forces. Since the unperturbed fluid is in equilibrium, the 
total radial force per unit mass F{r) at every point of the flow is zero, viz. 


F(r) 



r pr 


1 d 
p dr 


Pfi + 


Bl + Bt 


+ fl'r — 0 , 


( 1 ) 


where is the azimuthal flow velocity; and are the toroidal and axial components 
of the magnetic held in such units that the magnetic permeability is p = dvr; p and Pfi 
are the density and internal nonmagnetic pressure of the ideal huid, respectively; and Pr = 
—d^griw/dr is the radial force per unit mass due to the overall gravitational potential <hgrav 
Eq. (1) accounts for all the forces which act radially on each huid element in equilibrium; 
the centrifugal force, the tension of the magnetic held, the huid and magnetic-held pressure 
gradients, and the gravitational force. 

The change in free energy AE of the conhguration is dehned to be equal to the work 
done in performing the interchange between the positions ri and r2, viz. 


AE = - 


mF{r)dr — / mF{r)dr — / mPdV — / mPdV , 


' 1^2 


' 2^1 


' 1^2 


'2^1 


( 2 ) 


where the notation i ^ j with i,j = 1, 2 is used to indicate that the integral is taken during 
the displacement of huid element i from position i to position j] m = pV is the element mass; 
the change in volume dV = 0 in this method, as shown explicitly below; and P denotes the 
total internal pressure of each displaced huid element, i.e. 

Bl + Bl 

P = Pfl + Pmag = K p'' + , (3) 


where Pmag is the magnetic pressure, is a factor that depends on the specihc entropy of 
the huid, and 7 is the adiabatic index. K is constant for homoentropic huids and remains 
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the same during slow adiabatic interchanges, but varies in space if an entropy gradient is 
present in the fluid. 

Mass conservation in this method is expressed by the requirement that the total mass 
of the two fluid elements be the same before and after the interchange. Suppose that the 
interchange involves two fluid elements initially at radii ri and r 2 > ri with masses mi = piVi 
and m 2 = P 2 h 2 , respectively. After the interchange, the two masses are m'^ = p\V 2 and 
^2 = P' 2^1 radii r 2 and ri, respectively. Since the densities change adiabatically, i.e., 
p'l = Pi + ^Pad and P 2 = P 2 — ^Pad to linear order, then the conservation of total mass, 
mi + m 2 = m'l + m 2 , leads to the constraint (Ap — Apad)(y 2 — Vi) = 0, where Ap = p 2 — pi. 
Since, in general, Ap 7 ^ Ap^^, equal volumes must be adopted for adiabatic interchanges of 
individual fluid elements. 


It is now easy to understand the physics of the interchange method: During and after 
the interchange of two fluid elements, the internal pressure is not a driving force; it simply 
continues to adjust in trying to maintain the continuity of the fluid and the conserved quan¬ 
tities. Therefore, the pressure-gradient term in eq. (1) does not contribute to the integrals of 
eq. (2). On the other hand, the plus or minus signs of the purely radial forces are not changed 
by the interchange and these forces maintain the same unidirectional behavior. Only their 
magnitudes are changed when these terms are properly constrained by the applicable con¬ 
servation laws during the displacements from radius r* to Vj. Therefore, the radial forces will 
continue to drive the interchanged fluid elements either back toward their initial equilibrium 
positions ^restoring forces) or away from them (instability forces) and their combined con¬ 
tribution (in both magnitude and direction) to the free energy decides whether the system 
will evolve uphill or downhill in energy space. Of course, in the former case the system is 
bound to return back to the initial state of lower energy, while in the latter case it runs away 
from the initial equilibrium and toward nonequilibrium states of lower total energy. 

Since Vi = V 2 = V in the interchange method, we are then left with the following 
expression for the change in free energy written now per unit volume: 


AE 

'V 






dr. 


(4) 


As was mentioned above, a necessary and sufficient condition for stability to interchange of 
equal-volume elements is that 


AE 

'V 


> 0 . 


(5) 


In the incompressible case, the fluid density p is constant and the results of Paper I are 
obtained by simply evaluating the two integrals after the relevant conservation laws have 
been built into the force terms of each integrand. We now proceed to obtain the relevant 
stability criteria for compressible fluids in the two fundamentally different (due to differences 
in the conserved quantities) cases of interest (cf. Paper I): a Couette flow with a purely 
toroidal magnetic held, and a Couette how with both toroidal and axial held components. 


3. Purely Toroidal Magnetic Fields 

As was already strongly emphasized in Paper I, before we can calculate the above 
integrals, we must take into account all the physical quantities which are conserved during 
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the interchange. Clearly, total mass is conserved during the interchange (§ 2). Furthermore, 
because only a toroidal magnetic held B^{r) is assumed to be present, the specihc (i.e. 
per unit mass) angular momentum and azimuthal magnetic hux are also conserved during 
the interchange (see also Paper I).^ Finally, we assume that each huid element conserves its 
specihc entropy during interchanges. Note that, although the variations of the free energy are 
expressed per unit volume for convenience, the conserved quantities must be expressed per 
unit mass because individual element masses are not conserved during interchanges. Then, 
the above conservation laws are expressed mathematically by the following set of equations: 


mi + mj 

Li = [rv^]i 


K, 


.Pr\^ 

'Pfi 

cfi . 


const, 

const, 

const. 


const. 


( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


where and mj are the initial masses of the two elements; and the symbols L*, <Fi, and 
Ki refer to the conserved quantities in each huid element i (specihc angular momentum, 
specihc azimuthal magnetic hux, and specihc entropy, respectively). At this point, it should 
be noted that, owing to the various gradients present in the huid, the values of Lj, <hj, and 
Ki are generally different between any two huid elements located at two diherent radii r* 
and rj in the equilibrium how. 

Before we go on to obtain the stability criterion, we dehne what we mean by “adiabatic 
gradients” in huids threaded by a toroidal held only. An arbitrary perturbation in density 
dp leads to an adiabatic change in pressure dPad if the change occurs under the constraints 
imposed by the applicable conservation laws (iF=const and <h=const). In our case 


P = Kp^ + , 


( 10 ) 


and 


(W 

dr 


ad 


dr r 


( 11 ) 


The last term is written as B'^/r instead of rp^^"^ for convenience. The symbol represents 
the fast magnetosonic speed in the f direction: 


c 


2 

4> 


iPfi + Bl 
P 


( 12 ) 


The above dehnitions generalize the well-known hydrodynamic expression dPfi/dr\ad = 
c^dp/dr in which the adiabatic sound speed Cq in the huid is given by = '-yPfi/p. 


■^We note that, under the axisymmetric conditions assumed throughout § 3 and § 4, specific angular 
momentum is identical to circulation. For nonaxisymmetric perturbations, however, it is the circulation, 
and not the specific angular momentum, which dictates a conservation law. 



In a similar fashion, an arbitrary pertnrbation in pressure dP leads to an adiabatic 
change in density dpad if the change occurs under the constraints imposed by the applicable 
conservation laws (-h'=const and <h=const). Then we hnd that 


^ _ ^2 ^ I 

dr ^ dr r 


Combining now eqs. (11) and (13), we hnd that 


dP 

dr 


dP 

dr 





(13) 


(14) 


The magnetic held enters this equation implicitly through and the pressure terms. We 
shall see in § 4 that this behavior does not generalize in the presence of both and 
components, thus eq. (14) is not generally applicable to all magnetic-held conhgurations. 

We are now ready to calculate and reduce the integrals shown in eq. (4). Using the 
trapezoidal rule to evaluate each of the interchanges 1 —2 and 2 — 1 , we hnd that 


+ 


1 

2 

1 

2 


V 

^ + ’"2(Pl + ^Pad)^*hf} — {PrlPl + fi'r2(Pl + ^Pad)} 

^ I + {^2P2‘^2 + ^l(P2 — ^PadY^l} ~ {fl'r2P2 + 9rl{P2 “ ^Pad)} 


i+Ar) 

(—Ar). 
(15) 


Keeping terms up to second order in the expansions of eq. (15), eqs. (5) and (15) give 


AE 

'V 


pAr 


\AL‘^ 


— prA<l>^ 




Ap - Apad 
p 


>0 , 


(16) 


where A-quantities represent centered ”2 minus 1” diherences (e.g., AL?' = — and the 

remaining quantities represent averages over positions ri and r 2 > ri (e.g., m = (mi+ 7712 )/2). 
Using eqs. (1) and (14), eq. (16) becomes equivalent to the necessary and sufficient criterion 
for stability to axisymmetric perturbations 

1 dL^ _ d4>2 _ 1 UP _ BU up _ dP 

r^ dr ^ dr p'^c^ \dr r J \ dr dr 



with L = rv^ and <h = B^/{pr). 

Eq. (17) constitutes the hrst important result of this work. It is interesting that, apart 
from the new magnetoconvective term, eq. (17) has basically the same structure as the 
stability criterion derived for the much simpler incompressible case (see eq. [3.5] in Paper I 
where the constant p was normalized away). Of course, the above stability criterion also 
involves implicitly derivatives of the equilibrium density, a quantity that plays no role in the 
incompressible case. 
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The new magnetoconvective term in eq. (17) reveals the coupling between the mag¬ 
netic tension and the pressure differential that drives buoyancy: In the absence of held, 
this term simply introduces the hydrodynamic Brunt-Vaisala frequency Ng given by the 
expression (Tassoul 1978) 


IdPfi 

f 1 dp 

1 dPfA 

1 dPfi 1 

(dPfi 

dPfi 

) 

p dr 

Si- 

-8 

'yPfi dr ) 

dr ' 

\ dr 

dr 

ad/ 


(18) 


but in the presence of convective motions in astrophysical huids are actually strengthened 
by magnetic tension. This is because in gravitating huids, pressure gradients are intrinsically 
negative and B'^/r enters the second coefficient of the magnetoconvective term also with a 
negative sign. The result of magnetic tension then is to increase the magnitude of the last 
term of eq. (17) without ahecting its sign. Consequently, both the frequencies of stable, 
oscillatory waves and the growth rates of unstable magnetoconvective modes (if any) are 
increased by the B'^/r term accordingly. The behavior of the magnetic tension in this case 
is determined solely by the imposed conservation of specihc azimuthal magnetic hux at any 
given radial position r. 


Finally we note that the magnetoconvective term does not disappear in homoentropic 
huids with dK/dr=0 because total pressure depends implicitly on the magnetic held. For 
the present case with B^ held only, using eq. (10) and the conservation laws, the pressure 
diherential of eq. (14) can be written in the equivalent form 


dr 


dP 

dr 


dK 


ad 


_ 7 — - 1^2 

“IT + i'' IT 


(19) 


This form is used in the Appendix where we combine eqs. (17) and (19) with dK/dr=0 to 
conhrm Newcomb’s (1962) homoentropic stability criterion. 


4. Toroidal and Weak Axial Magnetic Fields 


The conservation laws of total mass and specihc entropy endure the introduction of an 
axial magnetic held Bz{r). On the other hand, the axial held eliminates the conservation 
laws of specihc angular momentum and azimuthal magnetic hux, and replaces them with the 
conservation of angular velocity and of axial current along held lines. A detailed description 
of this fundamental change was given in Paper I. In addition, the axial held introduces the 
conservation of the specihc axial magnetic hux. The relevant conservation laws of this case 
are expressed mathematically by the following set of equations: 


A 


+ rrij 

= const. 

(20) 

[’'Bd, 

= const. 

(21) 

'Vf 


= const. 

(22) 

- r 

i 


r BP 



(23) 



= const. 

. P . 

i 


\Pfi] 


= const. 

(24) 

.P\ 

i 
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where the new symbols /*, hi*, and Hi refer to the new conserved quantities in each fluid 
element i (axial current interior to the radius r,^ angular velocity, and specihc axial magnetic 
flux, respectively). 

Along the lines discussed in Paper I, we note that an axisymmetric perturbation in a 
fluid threaded by an axial held will generally induce curvature to the axial held and will 
therefore introduce curvature-dependent terms in the equations. The interchange method 
that we are implementing does not include such perturbations since 2 :-dependent terms are 
not accounted for during interchanges. On the other hand, if the axial magnetic held is 
sufficiently weak, its contribution to the energy change of the conhguration will be negligible 
compared to the other terms. Therefore, we conclude that the condition for stability which we 
are about to derive is necessary and sufficient only if the axial wavelength of the perturbation 
is inhnite, or in the limit Bz ^ 0. In general, this condition will be only sufficient for stability, 
since the curvature-dependent terms have a stabilizing inhuence (i.e., they make a positive 
contribution in the expression for AE). 

Having identihed the relevant conservation laws in this case, we can proceed just as in 
§ 3. First we dehne what we mean by “adiabatic gradients” in huids threaded by toroidal 
and axial held components. An arbitrary perturbation in density dp leads to an adiabatic 
change in pressure dPad if the change occurs under the constraints imposed by the applicable 
conservation laws (iF=const, /=const, and i?=const). In this case 


and 

dP 2 dp 

dr ^dr r 

The last term is written as H|/r instead of P/r^ for convenience. Notice the change of sign 
in this term relative to the corresponding term in eq. (11); it is caused by the introduction 
of axial current conservation in place of azimuthal magnetic hux. The symbol Cz represents 
the fast magnetosonic speed in the 2 : direction: 


(25) 

(26) 


.2 _ iPfi + Bl 

— 

p 

This is the characteristic speed in this case despite the presence of two held components. 
This fact alone tells us that the dynamics of the system is now driven by the magnetic 
coupling of huid elements in the 2 direction. 

In a similar fashion, an arbitrary perturbation in pressure dP leads to an adiabatic 
change in density dpad if the change occurs under the constraints imposed by the applicable 



^To be more precise, in order for the interchange to be performed without perturbing the background fluid 
and its magnetic fleld, screening currents must be set up at the radial surface boundaries of the displaced 
fluid elements. As these screening currents are moved along with the fluid elements, the total axial current 
interior to the radius r of each fluid element is indeed conserved during the displacement. 
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conservation laws (i^=const, /=const, and if=const). Then we find that 


dP 2 dp 

dr " dr r 


Combining now eqs. (26) and (28), we find that 


dr 


dr 





(28) 


(29) 


The two field components enter this eqnation implicitly throngh the pressnre terms. Bz is 
however the only field component that determines which magnetosonic speed is relevant to 
the stability of the given confignration. 

We are now ready to calcnlate and rednce the integrals shown in eq. (4). Using the 
trapezoidal rule to evaluate each of the interchanges 1 —> 2 and 2 —1, we find that 


AE 

'V 


1 

2 

1 

2 


+ (Pl + ^Pad)'/"2^l} + 

{P2^2^\ + {p2 — Apad)Tl^‘^'] + 



{PrlPl + 9r2i.Pl + ^Pad)} 
{9r2P2 + 9rliP2 ~ ^Pad)} 


i+Ar) 


(—Ar). 
(30) 


Keeping terms up to second order in the expansions of eq. (30), eqs. (5) and (30) give 


AE 

'V 


pAr 


rA^f 


AP 

prS 


+ Pr) 


Ap - Apad 

P 


>0 , 


(31) 


where A-quantities represent centered ”2 minus 1” differences (e.g., Afl^ = VL^ — Vl\) and the 
remaining quantities represent averages over positions ri and r 2 > ri (e.g., p = (pi + p 2 )/ 2 ). 
Using eqs. (1) and (29), eq. (31) becomes equivalent to the sufficient criterion for stability 
to axisymmetric perturbations 


^dD? _ ^ dP _ 1 UP BU up _dP \ ^ 

dr pr^ dr p^P, \dr r J \ dr dr ~ 


(32) 


with Q = v^/r and I = rB^. 

Eq. (32) constitutes the second important result of this work. Once again, apart from 
the new magnetoconvective term, eq. (32) has the same structure as the stability criterion 
derived for the much simpler incompressible case (see eq. [4.5] in Paper I where the constant 
p was normalized away). 

The new magnetoconvective term in eq. (32) reveals a change in coupling between the 
magnetic tension and the pressure differential when Bz is present (cf. eq. [17] in § 3 above): 
The magnetic-tension force now works to oppose buoyant motions and the development of 
convection in astrophysical fluids. This is because in the presence of gravity, pressure gradi¬ 
ents are intrinsically negative but B‘^/r enters the second coefficient of the magnetoconvective 
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term strictly with a positive sign. The difference has its origin in the presence of the Bz held 
that couples huid elements in the direction. When such elements attempt to move radially 
in opposite directions, the magnetic tension acts as a restoring force and tries to suppress 
these divergent motions. On the other hand, no such axial coupling exists between elements 
in the direction in the case of only held, where the behavior of the magnetic tension is 
determined solely by the conservation of specihc azimuthal magnetic hux (see § 3). 

Finally we note that the magnetoconvective term does not vanish in homoentropic 
huids with dK/dr=0 because total pressure depends implicitly on the two magnetic-held 
components. For the present case with B^ and B^ helds, using eq. (25) and the conservation 
laws, the pressure diherential of eq. (29) can be written in the equivalent form 


dP 

dr 



P’ 


dK 

dr 


1 

2 




(33) 


This form is used in the Appendix where we combine eqs. (32) and (33) with dK/dr=Q and 
I = 0, B^ = 0 to conhrm Newcomb’s (1962) homoentropic stability criterion. 


5. Static Magnetized Nonhomoentropic Atmospheres 


We are now in a position to discuss clearly the stability of static plane-parallel at¬ 
mospheres with and without magnetic helds. This is a topic of considerable contemporary 
interest due to its astrophysical applications to solar/stellar convection zones, stellar atmo¬ 
spheres, and sunspots (e.g., Thomas & Weiss 1992); and to accretion-disk magnetic helds 
(e.g., Frank, King, & Raine 1992; Duschl et ah 1994). Our results for compressible hows 
can be easily reduced to this special case by ignoring curvature-dependent terms (i.e., the 
centrifugal and the magnetic-tension forces); and by replacing the derivatives with respect 
to r with derivatives with respect to z, where Cartesian coordinates (x, y, z) are adopted 
throughout this section. In this new coordinate system, represents the vertical height in a 
plane-parallel atmosphere or convection zone while the gravitational force Qz = —d^g^av/dz 
points “downward” toward the negative ; 2 -axis. 


In the simplest case of a static, unmagnetized, plane-parallel atmosphere, the condition 
for stability is given by the convective term in eq. (17) or in eq. (32) with B^ = = 0, i.e., 

for stability to purely hydrodynamic convection 


dPfi f dPfi 
dz \ dz 


dP, 




dz 


ad 


> 0 


(34) 


where the equilibrium density gradient is given by the equation dPfi/dz = pgz- We note that 
the pressure gradient dPfi/dz < 0 since Qz < 0, and eq. (34) then leads to the Schwarzschild 
(1906) criterion (cf. Clayton 1983) 


dPfi ^ 'yPfi dp 

dz ~ p dz ^ 


(35) 


dPfi 

< 

dPfi 

dz 


dz 


( 36 ) 


or equivalently 
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The above analysis can be generalized now to inclnde the inflnence of a horizontal mag¬ 
netic held B = B{z)e:xi in the direction of the nnit vector Cx- In order to apply the interchange 
method correctly, we have to consider separately two types of vertical perturbations in an 
atmosphere: those with wavevectors k perpendicular to B (§ 5.1) and those with k parallel 
to B (§ 5.2). 


5.1. Pure—Interchange Instability: k T B 


We consider hrst the so-called “pure interchange” case in which the wavevector k T B 
(Tserkovnikov 1960; Newcomb 1961). In this limit, our analysis of § 3 is directly applicable. 
Replacing r hy z and by B^, dehning the conserved specihc hux by $ 3 , = B^/p, and 
ignoring curvature-dependent terms in the equations of § 3, the stability criterion becomes 
(cf. eq. [17]) 


where P = Pfi + p‘^^l/2 and 


with 


dz V dz dz 


dP 

dz 


-4 


ad 


2 _ iPfi + B, 


- ) >0 , 

(37) 

ad/ 


dp 
dz ’ 

(38) 

K 

(39) 


P 


For a nonhomoentropic 
the stability criterion 


atmosphere with < 0, then dP/dz < 0 and this result reduces to 


dP 

< 

dP 

dz 


dz 


(40) 


and using eqs. (38) and (39) along with the equilibrium expression dP/dz = pg^, we obtain 
the equivalent form (Tserkovnikov 1960) 


dp 

^ P^\9z\ 

dz 

- lP,i + Bl 


(41) 


This criterion for stability to pure interchanges is the direct generalization of the purely 
hydrodynamic Schwarzscliild criterion (eq. [35] or eq. [36]) and reduces to it in the limit of 
B.^0. 

Tserkovnikov’s (1960) result for k T B appears to suggest that the additional introduc¬ 
tion of a nonzero perturbation with k || B will not modify the sufficiency of this criterion 
since the restoring forces from held-line bending are naively expected to make a stabilizing 
contribution. This idea was sharply criticized by Newcomb (1961) who considered explicitly 
perturbations with k || B and obtained a different stability criterion which does not reduce 
to Tserkovnikov’s result in the limit of /c —0. 


This disagreement is simply a manifestation of structural instability between the two 
models. As was explained in Paper I, structural instability is caused by any physical process 
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that has the power to destroy/replace/modify one or more conservation laws. We show this 
in § 5.2 below, where we use the interchange method to examine the instability caused by 
perturbations with k || B. 


5.2. Parker Instability: k || B 

We consider now perturbations in the so-called “quasi-interchange” case (Newcomb 
1961). The interchange involves two horizontal thin hnite “strands” of fluid located initially 
along the magnetic held at heights Zi and Z 2 > Zi. We work in the limit of /c —0; hence, 
the derived criterion will only be sufficient for stability (see related discussion in § 4). 

The quasi-interchange case is more subtle than the pure-interchange case and its analy¬ 
sis proceeds along lines analogous to those discussed in § 4 above and in § 4.1 of Paper I: We 
allow for k ^ 0 perturbations by assuming that the huid and the magnetic held will remain 
unperturbed at a; = ±cx). The physical interpretation of this assumption is the following: 
During the displacement of huid element i from height Zi to height Zj (where i,j = 1,2), its 
mass is not conserved because it can how along the distorted held lines (i.e., in the direction 
of k) toward or away from the anchored (unperturbed) parts of the strand at a; = ±cx). 
Thus, just as in the i^^-held cases of § 4 above and in Paper I, the “horizontal inhnity 
X = ±cx)” serves in the present model as a “reservoir” through which mass is exchanged 
between perturbed huid elements. 

Under the assumption that the boundary conditions at a; = ±cx) do not change, the 
quantity that dictates a conservation law is now the Bernoulli integral in the direction of k 
(as found from the poloidal projection of the Euler equation; see, e.g., Lovelace et ah 1986): 

J = [ ^ + 4>grav = + 4>grav • (42) 

P iir=const 9” 


The above description is the essence of the so-called Parker (1966, 1975, 1979) instability 
(see also Shibata et ah 1989) for which we are now equipped to apply the interchange 
method. The relevant adiabatic density gradient in the quasi-interchange case is determined 
from eq. (42) and the new conservation law dJ/dz = 0, viz. 


dz ad ~ iPfi 


(43) 


The magnetic held does not enter the Bernoulli integral because k || B perturbations 
evolve only along held lines where the magnetic held does not exert tension. This again 
hints that we should interchange two huid elements of diherent masses rrii ^ rrij and equal 
volumes Vi = Vj. By considering the interchange of the contents of two equal volumes, the 
conservation of specihc magnetic hux is expressed by the requirement that the held strength 
remain the same per unit mass; and, as was shown in § 2, the conservation of the total mass 
contained in the two volumes is guaranteed. 

Ignoring then curvature-dependent terms and keeping terms up to second order in the 
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calculation of the gravity integrals of eq. 


(4), we obtain from eq. (5) the condition 


+ 9z2{,Pi + ^Pad)] 

+ 9zi{p2 - Apad)] (-A^:) > 0, (44) 

or 

A£/ 

= Az [gz{Ap - Apad)] > 0 , (45) 

where A^ = Z 2 — zi and pz = {pzi + gz 2 )/‘^- Since A^ > 0, eq. (45) is equivalent to the 
sufficient criterion for stability 


AE 

'V 


- - [gziPi 


7^ [gz2P2 


(‘k.'k )>o 

\dz dz 


(46) 


Recalling that Pz < 0 and using eq. (43), the stability criterion can be written in the general 
form (Newcomb 1961) 


^ > P^\9z 
dz ~ 'jPfi 


(47) 


This inequality expresses the sufficient condition for suppression of the Parker instability in 
the limit of /c —0 perturbations. The stability criterion does not depend on the magnetic 
field and is thus different from eq. (41). Thus, it is now clear from our analysis that structural 
instability (the change in the conservation laws) explains the difference of results in the cases 
with k T B and k || B. 


6. Summary and Discussion 

We have derived general forms of the axisymmetric stability criteria in compressible, 
rotating, magnetized, nonhomoentropic Couette flows using an energy variational principle, 
the interchange or Chandrasekhar’s (1960) method. The results given in § 3 and § 4 gen¬ 
eralize previous limited expressions obtained for incompressible flows (Chandrasekhar 1981; 
Paper I), for magnetostatic models (Tserkovnikov 1960), and for suppression of magnetorota- 
tional instabilities in compressible homoentropic flows (Newcomb 1962; Rogers & Sonnerup 
1986). To demonstrate the generality of our results and to assert the power of the inter¬ 
change method, all these previous results are rederived as special cases of our two general 
stability criteria (eqs. [17] and [32]) and are placed in context in an Appendix to the paper. 

We should point out that our resuls have been obtained by considering adiabatic per¬ 
turbations (thus also frozen-in magnetic fields). Our strongest and most general stability 
criterion (32) differs from those obtained in magnetized models with finite conductivity and 
equal initial temperature along field lines, such as the criteria of Balbus (2001). In our 
treatment in § 4, perturbations are constrained to conserve angular velocity, axial current, 
and specific entropy; it is precisely the gradients of these physical quantities that appear in 
the individual terms that make up the stability criterion (32). On the other hand, the linear 
analysis of Balbus (2001) is based on isothermality along field lines and the Boussinesq 
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approximation, conditions that impose the conservation of temperature between equally hot 
elements coupled by the same held line. As a result then, temperature replaces specihc 
entropy in that stability criterion. A potential problem in this approach is that coupled 
perturbed huid elements that are maintained at the same temperature by the new conser¬ 
vation law are, at the same time, not allowed by the Boussinesq assumption to change their 
pressures. Nonetheless, the structural instability discussed at length here and in Paper I is 
ever present in this model too, owing again to the change in conservation laws effected by 
the axial held; as Balbus also points out, the hydrodynamic limit (the Hqiland criteria; see 
Tassoul 1978) can never be recovered even if the magnetic held is allowed to vanish. 

We believe that our results are relevant in several areas of contemporary astrophysical 
research. We outline below three areas along with brief discussions of the pertinent issues: 

(a) Magnetized stellar convection zones: The interchange method is directly applica¬ 
ble to static, magnetized, plane-parallel atmospheres and stellar convection zones (see § 5), 
where it recovers well-known criteria for suppression of convective motions and magnetic 
buoyancy (Schwarzschild 1906; Tserkovnikov 1960; Newcomb 1961; Parker 1966, 1975, 1979). 
In addition, this method captures easily and explains clearly the diherence between previ¬ 
ously obtained results (Tserkovnikov 1960; Newcomb 1961) by identifying the apparent dis¬ 
agreement between magnetoconvective stability criteria (eq. [41] and eq. [47]) with a typical 
case of structural instability. 

(b) Magnetized low-luminosity accretion flows: We believe that our work also sheds 
some light to the recent controversy concerning the radiatively inefficient accretion hows and 
the direction of angular momentum transport in associated models (Narayan et al. 2000, 
2002; Balbus & Hawley 2002). Our analysis of § 4 delineates the role of rotation, magnetic 
helds, and compressibility to the stability of such hows. As emphasized here and in Paper I, 
how stability as well as angular momentum and magnetic hux transport are governed by the 
demand that certain conservation laws be respected by perturbed huid elements. Specihcally, 
for elements coupled by the same poloidal held line: 

1. The constant O = L/r^ law leads to outward angular momentum transport—as some 
huid elements move outward their r increases, thus their L must also increase in order 
to maintain the same H. 

2 . At the same time, the constant I = rB^ law leads to inward azimuthal magnetic hux 
transport—as some huid elements move inward their r decreases, thus their must 
increase in order to maintain the same I. 


The sense of transport of these quantities is absolutely hxed by the mere presence 
of an axial magnetic held component, the component that is also directly responsible for 
the existence of magnetorotational instability. The magnetized how model of § 4 makes it 
quite clear that the non-magnetic stability criteria and the sense of convective transport of 
associated quantities cannot be smoothly recovered, even in the limit of 0. It is thus 

hard to see how angular momentum could be transported inward in the presence of a frozen- 
in axial held. This however might be the case if the held could “slip” through the plasma 
(in which case the ideal MHD conditions do not apply), a situation that would immediately 
also invalidate the conservation laws (21) and (22). 

In relation to the role of the B^p component, we re-iterate here that, in the presence of 
an axial held, this component can reduce signihcantly the ehects of bouyancy, as shown by 
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eq. (32) above. It is unfortunate that the component has been ignored in many previous 
studies with no good justihcation, other than that its absence makes the linear-stability 
analyses manageable. We believe that such simplihcation of the dynamical equations is not 
justified in accretion disks for the following reasons: (i) Differential rotation will twist around 
any initial weak radial field and will continue to amplify the resulting toroidal held, (ii) The 
toroidal held produced in this way will be transported inward, just as the angular momentum 
is transported outward. It is therefore certain that azimuthal magnetic hux will pile up in 
the inner regions of magnetized disks where it will dominate the dynamics of these regions. 
We discuss this point in more detail in part (c) below. 

While the instability criteria that we have derived are ehectively local, it is interesting 
to speculate about possible applications of our results to global hows, in particular those of 
ADAF-type, in which all sound speeds scale as c* oc and the temperature scales as 

T oc r~^ (Ichimaru 1977; Narayan & Yi 1994). A similar issue was raised in Paper I, where 
we speculated about whether nonlinear simulations might lead to hows with constant D and 
oc 1/r that hover around the incompressible Ai7 ~ 0 stability limit. Now, it is clear from 
the outset that D cannot be constant in an ADAF. However, the formulation of our problem 
still allows us to examine the possibility of having a how with Ai7 ~ 0 on all scales, i.e. a 
how that operates at the edge of instability at all relevant scales. Then the condition AE ~ 0 
requires that all terms in eqs. (17), (32) be of the same order of magnitude; replacing the 
d/dr derivatives by 1/r in either equation, we hnd that 







(48) 


over all scales, i.e., that the magnetic held should be in rough equipartition with the rotational 
kinetic energy! This condition is independent of the specihc radial scalings obeyed by D(r), 
B^{r), and p(r). 

It is interesting to note that the global condition (48) obtains only in the compressible 
case, since in an incompressible how the density drops out of the stability criterion (cf. 
eq. [4.6] of Paper I). It is also interesting that such equipartition was actually assumed in the 
original ADAF formulation for the explicit purpose of reducing the adiabatic index of the 
how to below 5/3. Finally, we point out that although eq. (48) does not uniquely determine 
the scalings of D(r), B^{r), and p(r) of the how, it still imposes an important constraint on 
the strength of the toroidal held for a given accretion geometry. 

(c) Magnetized jet-like outflows: Our results are useful to the analysis of the vertical 
structure of fully ionized magnetized accretion disks. The stability criterion (41) may resolve 
the issue of how much azimuthal magnetic hux can be stored in the interior of an accretion 
disk before it begins to escape buoyantly to the surface and, from there on, to the surrounding 
region. The condition for stability against buoyant escape of a whole azimuthal magnetic-hux 
tube (eq. [41]) is weaker than the condition against the “nonaxisymmetric” Parker instability 
(eq. [47]) which requires steeper density gradients to stabilize k || B perturbations. However, 
if we consider axisymmetric perturbations in an axisymmetric accretion disk with only a 
toroidal-held component, then only eq. (41) is relevant and the criterion for stability against 
vertical magnetic buoyancy in an homoentropic disk can be written in the form (cf. eq. [19] 
with dK/dr = 0 and without curvature-dependent terms; and eq. [37] with dP/dz < 0) 

dz dz V p^ / 


> 0 . 


(49) 
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A similar condition for suppression of radial magnetoconvective motions, namely 


dH^ 

dr 


d 

dr 



>0 , 


(50) 


is obtained from eqs. (32) and (33) in the special case of a nonrotating (hi = 0) homoentropic 
{dK/dr = 0) disk with a dynamically important, purely axial (J = 0, = 0) magnetic 

field Bz, and with dP/dr < 0 (see also Anzer [1969] and Lubow & Spruit [1995]). Both 
conditions are expressions of the Kruskal-Schwarschild (1954) criterion for stability of an 
atmosphere supported against gravity only by fluid-pressure and “horizontal” magnetic- 
pressure gradients. 

We heretofore focus on the signihcance of eq. (49) for jets and accretion disks in active 
galactic nuclei. It is generaly believed that the toroidal magnetic held, which is naturally 
produced by differential winding of any weak radial magnetic-held component, is always 
limited by buoyant escape to remain much below equipartition levels (see, e.g., D’Silva & 
Chakrabarti 1994). On the other hand, numerical simulations (Shibata, Tajima, & Mat- 
sumoto 1990) indicate that under some (unspecihed) conditions, buoyant escape appears to 
be suppressed and does not prevent the toroidal magnetic held from growing to equipartition 
values (see also Galeev, Rosner, & Vaiana 1979). In the same context, Contopoulos (1995) 
has considered the role of the toroidal magnetic held in generating collimated energetic bipo¬ 
lar jet outhows in accretion disks. The issue of the initial buildup of the toroidal held was 
however left as an open question. 

We now suggest the following scenario: Even if the held is initially conhned to the 
interior of the disk, as inhow of matter and diherential winding of the radial magnetic held 
continue, inequality (49) can be reversed. Then, magnetoconvective instability will set in 
and azimuthal magnetic hux will leak out of the disk, in an attempt to return to a stable 
vertical stratihcation with d^'^/dz > 0. Now, one has to consider the relevant timescale 
for radial inhow of matter and azimuthal magnetic hux; If the growth of in the disk’s 
interior proceeds dynamically (because of an implosive magnetorotational instability), then 
the buoyant escape will take place explosively—^just as during the release phase of a typical 
laboratory plasma gun—and will presumably result in the discharge of highly energetic, 
axially moving “blobs” of disk material. This mechanism may account for the relativistically 
moving blobs observed by VLBI in jet-like outhows in active galactic nuclei (e.g., Hughes 
1991; Burgarella, Livio, & O’Dea 1994). 

Condition (49) for suppression of vertical magnetoconvective motions is everywhere 
satished if the huid density decreases sufficiently fast (e.g. exponentially) with height from 
the mid-plane of an accretion disk. In contrast, if the density does not decrease sufficiently 
fast with height, then inequality (49) may be permanently reversed at some height. Such a 
reversal will lead to slow leakage of the inhowing azimuthal magnetic hux. This continuous 
leakage in the direction may prevent altogether the dynamical buildup of the azimuthal 
hux which is expected to produce by discharge a relativistic jet outhow (see the discussion 
above and in Contopoulos 1995). This tentative conclusion deserves further investigation in 
relation to the dichotomy between radio-loud and radio-quiet quasars. 


This work was supported in part by a Chandra Cuest Observer grant. 



A. Special Forms of the Stability Criteria 


A.l. Comparisons with the Results of Tserkovnikov (1960) 


(a) In order to compare our stability criterion for an homoentropic magnetostatic cylindrical 
equilibrium threaded by a purely toroidal magnetic field to the corresponding expression 
found by Tserkovnikov (1960) and confirmed by Newcomb (1962), we combine: 

1. Eq. (1) with ^ 0 = 0 , (7^=0, and - 02 = 0 ; 

2. Eq. (17) with L=0; 

3. Eq. (19) with dK/dr=0] 


and we obtain the form 


dr 


< 0 . 


(Al) 


After straightforward algebraic reductions, this stability criterion can be cast in the equiva¬ 
lent form 


d\n.Bs ^ iPfi - Bi 


< 


d\nr 'jPfi + B^ 


2 ’ 


(.42) 


that was derived by Tserkovnikov (1960). 

We see now how the complexity in Tserkovnikov’s result obscures the basic physics 
of the magnetostatic problem. Eq. (A2) is just a complicated relation between B^ and 
Pfi with no obvious physical meaning. On the other hand, the equivalent form eq. (Al) 
reveals that, despite the homoentropic assumption, unstable convective modes can still exist 
if the magnetic field establishes a superadiabatic pressure gradient. This result is valid 
in the incompressible case as well (Chandrasekhar 1981; eq. [3.7] in Paper I). Its general 
applicability should have been expected because the second-order pressure differential in 
eq. (19) depends also on the magnetic flux gradient (and not on the entropy gradient alone). 

(b) Along similar lines, it is expected that the stability criterion for a magnetostatic equilib¬ 
rium threaded by both toroidal and axial field is expressed by a simple condition involving 
the radial variation of the axial current. This result has been previously derived in the in¬ 
compressible case (Chandrasekhar 1981; eq. [4.7] in Paper I) but it has not been generalized 
for compressible models until now. We combine: 

1. Eq. (1) with V 0 =O and ( 7 r= 0 ; 

2. Eq. (32) with f2=0; 


and we obtain the form 


dP 

dr 


< 0 


(43) 


Note that eq. (33) plays no role here in interpreting the result because the reductions = 0 
and = 0 in eq. (1) imply that dP/dr + B'^/r = 0 and the last term of eq. (32) is 
then eliminated from the stability criterion. Therefore eq. (A3) is generally valid for any 
distribution of field Bz{r) and axial flux H{r) = Bz{r)/p{r) in magnetostatic models. 
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A.2. Comparisons with the Results of Newcomb (1962) 

The stability criteria obtained by Newcomb (1962) were also discussed in the Appendix 
of Paper I where we reported erroneously that Newcomb’s derivations miss an additional 
magnetoconvective term. This statement is not correct—it turns out that the slightest error 
in the application of the interchange method generates additional spurious terms. However, 
the mathematical reductions reported in the Appendix of Paper I were performed only on 
Newcomb’s equations and they are correct. 

Using the general stability criteria reported in this paper, we can reproduce now New¬ 
comb’s results precisely: 

(a) In order to compare our stability criterion for an homoentropic fluid threaded by a 
purely toroidal magnetic held to the corresponding expression obtained by Newcomb (1962), 
we combine: 


1. Eq. (1) with gr=0 and 5^=0; 

2. Eq. (17) with $ = B^/{pr-); 

3. Eq. (19) with dK/dr=0 and <h = B^/{pr); 


and we obtain the form 
1 d 

r3 


- 


2 dBf 
r dr 


+ 


2Bl 


> 


r2 ' pci 


pB^ - 2r^Bl 


m 


which is identical to eq. (8.11) in Newcomb (1962). We see again how the complexity in 
Newcomb’s result obscures the basic physics of the problem. Eq. (A4) is just a complicated 
relation between all the relevant quantities with no obvious physical meaning. On the other 
hand, the equivalent form of the stability criterion (17) along with the deflnition (19) can 
be physically interpreted with ease as the interplay between gradients in the three conserved 
quantities (speciflc angular momentum, azimuthal flux, and entropy), as was done in § 3 
above. 

(b) In order to compare our stability criterion for an homoentropic fluid threaded by a purely 
axial weak magnetic held to the corresponding expression obtained by Newcomb (1962), we 
combine: 


1. Eq. (1) with gr=0 and B^=0; 

2. Eq. (32) with 1=0 and B^=0] 

3. Eq. (33) with dK/dr=0, 1=0, and H = B^/p', 


and we obtain the form 


A 

dr 


{pL^) > 


ApB 


+ 


pel 



(A5) 


r 
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which is identical to eq. (10.4) in Newcomb (1962) if we let Hz —0 and Cz —> Co (see 
eq. [27]). Once again, eq. (A5) is just a complicated relation between all the relevant quan¬ 
tities with no obvious physical meaning. On the other hand, the equivalent form of the 
stability criterion (32) along with the definition (33) can be physically interpreted with ease 
as the interplay between gradients in the three conserved quantities (angular velocity, axial 
current, and specific entropy), as was done in § 4 above. 


A.3. Comparisons with the Results of Rogers &; Sonnerup (1986) 

Rogers & Sonnerup (1986) have obtained a form of the stability criterion in the case of 
an homoentropic fluid threaded by toroidal field only, and they pointed out that their result 
agrees with the condition found by Newcomb (1962), eq. (A4) above. Although technically 
correct, their result is subject to two fundamental assumptions: (a) The fluid is uniformly 
rotating, while Newcomb’s result is valid for a differentially rotating fluid, (b) The fluid 
is gravitating, while Newcomb’s result was derived for a nongravitating fluid. Therefore, 
according to (a), the criterion of Rogers & Sonnerup is only a special case of Newcomb’s 
result, but according to (b), it is more general than Newcomb’s result. Obviously, we must 
exercise caution when we reduce our result to the form of Rogers & Sonnerup (1986), and 
we should not use eq. (A4) because it was derived under the assumption that gr=Q- 

In order to compare our stability criterion for an homoentropic, uniformly rotating, 
gravitating fluid threaded by a purely toroidal magnetic field to the corresponding expression 
obtained by Rogers & Sonnerup (1986), we combine: 


1. The assumption that = const; 

2. Eq. (1) with = fir, Hz=0, and Qr = —d^grav/dr; 

3. Eq. (17) with L = flr^ and <l> = B^/{pr); 

4. Eq. (19) with dK/dr=0 and <l> = B^/{pr)-, 

5. The definitions of Rogers & Sonnerup (1986): Pp = Kp'^, Pmag = P = Pfi + 

Pmagi &Ild T = ^grav fl r /2, 


and we obtain the form 

2 / 2Pmag _ dPmag \ _ 1 / ‘^Pmag _ (^\ ^ ^^^2 > g 

r \ r dr J iPfi + ‘^Pmag \ r dr J dr dr ~ 


(^ 6 ) 


which is identical to eq. (15) in Rogers & Sonnerup (1986). We note, once again, that 
fl = const in this derivation, hence dT/dr = — [pr + fl^r) in eq. (A 6 ). 

In this case too, the above result contains substantially less physics than our stability 
criterion (17)—no gradients in the angular velocity or in the specific entropy—yet it is hard to 
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deduce the elegant physical interpretation afforded by our equivalent form (eqs. [17] and [19] 
with n = const and dK/dr=0): 


dr 2c^ \dr r J dr 


(717) 


Here the hrst term is the epicyclic frequency due to solid-body rotation and has a stabiliz¬ 
ing influence. The second term depends on the azimuthal flux gradient. The last term is 
convective and is amplified by the toroidal held, but its sign depends only on the sign of 
the gradient d^"^/dr. This is because dP/dr < 0 and thus the coefficient {dP/dr — B‘^/r) 
is always negative in cases of astrophysical interest. Consequently, the pressure differential 
is determined to be superadiabatic or subadiabatic solely by the radial distribution of <h^, 
given that the huid is homoentropic (eq. [19] with dK/dr = 0). 
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